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1.  Indroduction.  Let  be  independent  populations  represent¬ 

ing  k  experimental  treatments  and  let  ttq  be  the  control  treatment.  Let 
f(x,6.j)  denote  the  density  of  ,  i=0,l,...,k.  Any  population  iri  is  said 
to  be  superior  to  the  control  if  0 .  >  6q,  and  inferior  otherwise.  While 
00  is  not  known,  we  have,  based  on  past  experience,  a  fair  idea  of  it  so  as 
to  assume  that  eQ  <_  0g,  a  known  quantity.  Following  the  earlier  setup  of 
Gupta,  Huang  and  Nagel  [1]  and  Huang  and  Panchapakesan  [3],  who  have  stud¬ 
ied  locally  optimal  rules  based  on  ranks  for  selecting  the  best  population, 
we  assume  that  the  functional  form  of  f(x,e)  is  known  but  for  the  value  of 
the  parameter.  We  seek  a  procedure  based  on  ranks  in  view  of  the  usual 
considerations  of  robustness  against  possible  deviations  from  the  model. 

We  are  interested  in  selecting  a  subset  (possibly  empty)  of  the  k  experi¬ 
mental  treatments  consisting  of  those  that  are  superior  to  the  control. 

Let  X^,  j=l,...,n,  be  independent  observations  from  ,  i=0,l,...,k. 
Let  R^  denote  the  rank  of  X^j  in  the  pooled  sample  of  N  =  (k+l)n  obser¬ 
vations.  The  smallest  observation  has  rank  1  and  the  largest  rank  N.  Let 
X1  1  *2  -  -  XN  ^enote  tbe  ordered  observations.  A  rank  configuration 

is  an  N-tuple  A  =  (Ap...,A^),  A^  €  (1 ,2,. . .  ,kl,  where  A..  =  j  means  that  the 
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ith  smallest  observation  in  the  pooled  sample  comes  from  ir..  Let  c  -  (a) 

J 

denote  the  set  of  all  rank  configurations  for  fixed  k  and  n.  For  fixed  a, 
let  X  =  {x€2,|ax=a},  where  %  =  {x:  x=(x-j ,. . .  ,x^)}  and  Ax  denotes  the 
rank  configuration  of  x  *  (x-| ,. . .  .x^).  A  decision  rule  6  based  on  the  ob¬ 
served  rank  configuration  A  is  a  k- tuple  defined  by  6  =  <s(a) 

=  { 6 1  (a)  , . . .  ,  6^(a ) } ,  where  6^(a)  is  the  (conditional)  probability  that 
is  selected  as  a  superior  population. 

Let  o  =  (0q,6j  .  ,9^)  and  a  =  {0|0q1O*}.  Define 

nQ  =  ^|0i=9oi0o’i=1**--’k}  and  Q*0  =  {0l0j=eo< 

We  are  interested  in  the  class  of  rules  6  satisfying 
(1.1)  P0{n.j  is  selected|e  e Uq}  <.  y  for  i=l,...,k. 

In  this  class,  we  seek  a  locally  optimal  rule  in  the  sense  that  it  maximizes 


(1.2) 


k 

J  -g-^7  Pgliij  is  selected! e 


6i=0O 


Let  P  (a)  denote  the  probability  of  realizing  the  rank  configuration  A. 
Then  (1.1)  can  be  written  as 


0.3) 


l  6i (a)Pq  (a)  <  y  for  i=l ,. . .  ,k, 


^0 


where  eQ  =  (0q, . . . ,eQ)  € Hq  and  the  expression  (1.2)  is  equal  to 


(1.4) 


l  1  6i (A)p  (i )(A) 

i=l  30  i  c  1  on' 


where  e 


(i) 


denotes 


er8o 


a  point  in  n*Q.  The  condition  (1.3)  corresponds  to  controlling  error  prob¬ 
abilities  and  the  optimality  condition  in  (1.4)  reflects  the  sensitivity  of 
the  rule  when  all  but  one  population  are  not  distinctly  superior  (ej=0g,j^) 
and  the  remaining  one  is  in  a  neighborhood  of  the  others  but  distinctly 


superior  (o^  >  eg). 
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2.  Derivation  of  a  locally  optimal  rule.  We  assume  that  the  density 
f(x,e)  satisfies  the  following  set  of  regularity  conditions:  (i)  f(x,e)  is 
absolutely  continuous  In  e  for  almost  every  x,  (ii)  f(x,e)  Is  continuously 
differentiable  with  respect  to  e  for  almost  every  x,  and  (iii)  f(x,e) 

=  —  f(x,e)  is  integrable. 

Now,  the  probability  P „(a)  of  realizing  the  rank  configuration  a  under 
0  €  ft  is  by 


(2.1) 


2  N 


P0(A)  =  /  /  ...  /  n  f(xi ,e  )dxr..dxN  . 

—  —  OO  —  CO  -.00  i=l  ■) 


We  note  that  P.  (a)  is  independent  of  the  common  value  en  of  the  parameters 
2q  u 


and  is  equal  to  1/N!. 
(2.2) 


Thus,  the  condition  (1.3)  becomes 

jjr  1  A j  (A)  ±y  for  . k- 


For 


it  can  be  easily  seen  that 


(2*3)  ae7  P0(i)(A) 


51  Ai(A,e*),  say. 


Thus  we  wait  to  derive  a  rule  6  which  satisfies  (2.2)  and  which,  among  all 
rules  that  satisfy  (2.2),  maximizes 
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which  in  general  depends  on  og.  However,  it  is  independent  of  eg  if  it 
is  a  location  or  scale  parameter. 


3.  A  special  case.  One  can  specialize  the  rule  6Q  given  by  (2.5) 
to  specific  densities  f(x,e).  An  Important  special  case  arises  when  f(x,B) 
is  the  logistic  density  f(x,e)  =  e"^x‘e V[l+e“^x"eh2*  <  x  <  «  , 

<  e  <  ~  .  In  this  case,  4>(u,f,e)  =  2u  -  1  which  leads  to  equally  spac¬ 
ed  scores  and  is  of  the  form  =  a  +  ib,  where  b  >  0.  Consequently, 
the  rule  is  given  by 


(3.1) 


«°(a)  • 


1 

P 

l  0 


I  Rii  -  C/N! 
j=1 


where  0  <  p  <  1  and  c  are  determined  by 


(3-Z)  Pog  { jl,  Vc/N!lt‘'|,e5{j1  Vc/Nl)^' 

The  values  of  p  and  c  can  be  obtained  from  tables  for  Wilcoxon  two-sample 
rank-sum  statistic. 


4.  Some  remarks.  Nagel  [4]  defined  just  rules  for  selecting  the  best 
population.  This  concept  can  be  applied  also  to  the  problem  of  selecting 
populations  that  are  better  than  control.  In  our  setup,  it  means  that  the 
probability  of  selecting  is  nondecreasing  if  all  the  observations  from 
Tfj  are  increased  and  the  observations  from  all  other  populations  are  decreas 
ed.  The  rule  6°  defined  by  (2.5)  is  just  if  is  nondecreasing  in  i.  In 
the  case  of  location  parameters,  this  monotonicity  of  Bi  is  equivalent  to 
saying  that  f(x)  is  strongly  unimodel,  i.e.,  -  log  f(x)  is  convex  (see  [2], 
p.20).  In  the  special  case  of  logistic  densities,  the  rule  6°  given  by 
(3.1 )  is  just. 
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Though  eQ  is  not  known,  we  have  assumed  that  an  upper  bound  eg  is 
known.  If  eQ  is  known,  then  in  stating  the  optimality  requirement,  eg 
Is  replaced  by  oQ. 
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